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ABSTRACT 


Differential equations and boundary conditions for a rotor blade in 
forward flight, with subsonic or transonic tip Mach number, are derived. 

A variety of limiting flow regimes determined by different limits involv- 
ing blade thickness ratio; aspect ratio, advance ratio and maximum tip 
Mach number is discussed. The transonic problem is discussed in some 
detail, and in particular the conditions that make this problem quasi- steady 
or essentially unsteady are determined. Asymptotic forms of equations and 
boiindary conditions that are valid in an appropriately scaled region of the 
tip and an azimuthal sector on the advancing side are derived. These equa- 
tions are then put in a form that is valid from the blade tip inboard through 
the strip theory region. 
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SYMBOLS 


a = - .(^2 X r + V) 
c 

f = f(x. y) 


F=F(x, y, z, t)=0 
g = g(x, y) 
i, j 

i 

M = QR/c 
m = M(1 + u) 

P 

r 

R 


S = 


t 

T 

V 


X. y, z 


y 

6 

e. 

u 

T 


free stream velocity relative to blade fixed coordinates 
sound speed 

non-dimensional blade geometry in blade attached co- 
ordinate system 

blade geometry relative to an inertial coordinate sys- 
tern 

dimensional blade geometry in blade attached coordinate 
system 

unit vectors along x and y axes 
blade chord 

reference Mach number 

tip Mach number on advancing side 

static pressure 

position vector 

blade radius 

distance in chords measured toward center of rotation 
from blade tip 

span coordinate for transonic problem 
time 

scaled time variable for transonic problem 
forward velocity of rotor 

coordinates along the blade chord, span and normal 
to the plane of rotation 

ratio of specific heats 

scale parameter for time in transonic problem 
inverse of blade aspect ratio 
advance ratio 
thickness ratio 
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ip 


velocity potential 
reduced frequency 
rotational velocity 


w 

Q 


A prime refers to a fixed (inertial) coordinate system. 
Subscript oo refers to conditions in undisturbed air. 

A bar denotes a vector. 
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UNSTEADY SUBSONIC AND TRANSONIC POTENTIAL. FLOW 
OVER HELICOPTER ROTOR BLADES 

I. INTRODUCTION 

Helicopter rotor blades that could support a nearly shock free tran- 
sonic flow in the advancing blade tip region would improve helicopter aero- 
dynamic performance and would reduce noise produced by the blade. The - 
design of better high speed rotor blades has so far been guided mostly by 
qualitative consideration of how blade geometry influences the flow in the 
complex aerodynamic environment of the tip region. While such considera- 
tions are' undoubtedly useful, they necessarily involve some guessing of how 
blade geometry influences flow behavior, particularly at high speeds. 
Verification of a particular designmust rely on expensive (and time con- 
suming) wind tvuinel and flight testing. It is therefore clear that blade de- 
sign (as in any other design problem) would benefit from analytical and 
numerical methods that would predict the details of a transonic flow field 
in the tip region for a given blade geometry and specified flight parame- 
ters. Calculation of the nonlinear, three-dimensional and possibly un- 
steady flow near a rotor tip is, of course, only a first step toward better 
design. The problem here is that while there have been recent advances 
in the analytical and numerical methods that will generate shock-free two- 
dimensional airfoil sections, there are no such methods available for 
three-dimensional flows. There is also ho indication yet that a shock- 
free two-dimensional section would show better performance in a three- 
dimensional environment than any other standard airfoil section, although 

it may be that the two-dimensional designs do have advantages over stan- 

\ 

dard sections in a three-dimensional application. 
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Preliminary calculations of transonic pressure distributions on 
rotor blades were made by Caradonna and Isom (1). They considered 
what is probably the simplest case, a non-lifting, hovering rectangular 
blade with 6% thick biconvex sections operating in the lower transonic 
regime. The calculations show how direct numerical methods suitable for 
two-dimensional problems can be extended in a fairly straightforward way 
to the three-dimensional problem of steady transonic flow near a rotor 
tip. 

The analysis of reference (1) is restricted to a rotor with zero ad- 
vance ratio, so that the flow is steady in a coordinate system attached to 
the blade. In this report we shall consider the more complicated problem 
of a rotor with non- zero advance ratio. The primary objective is to ob- 
tain suitably scaled differential equations that describe transonic flow over 
the advancing rotor tip, and to show how possible flow regimes depend on 
various limiting cases of scale parameters. In particular, it is important 
to isolate what might be called the "worst" case, which is a three-dimen- 
sional, nonlinear and \msteady (as opposed to quasi- steady) flow regime on 
the advancing side of the blade. 

The hover problem of reference (1) is characterized by three para- 
meters: blade thickness ratio t, aspect ratio AR = e ^ and tip Mach 

number M = SIR/ c (S2 is angular velocity, R is blade radius and c the 
speed of sound in undisturbed air). It was shown in reference ( 1) that the 
hover problem is in some ways similar to the problem of flow over a fixed 
wing, and the analysis was gviided by this (sometimes rough) analogy. The 
advancing rotor problem involves one additional parameter, the advance 
ratio lu , as well as the time variable t. The reference Mach number is 
then the maximum free- stream Mach number on the advancing side. 
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JH = M(1 + ^). Our objective is to determine how T, e, JH and /u combine 
as coefficients in the governing differential equations and boundary condi- 
tions, and then to consider asymptotic forms of these equations, with t 
and e small, TH either subsonic or transonic, and u either order one or 
small. The asymptotic results always simplify \anwieldly equations, and, 
more importantly, they serve to put equations into a form that is subject 
to useful physical interpretation. It will be noted, however, that pursmt 
of a variety of limiting cases may produce distinctions that are rather fine, 
and in some instances too fine. For example, we consider three limiting 
values of advance ratio* a high advance ratio problem where ^ is of the 

2/3 

order of 0. 4 and fixed; a case for which ju = 0 (t )j and one for which 
/n = 0 (t), with T -♦ 0 in each case. Each of these limiting cases results in a 
distinct set of governing equations and a distinguished flow physics. How- 
ever, when the blade thickness ratio t is about 0. 1, then /u = 0. 4, or 

2/3 

H = 0(t ' ) or u = 0(t) are not widely different values for advance ratio. 

Nevertheless, the usual scaling procedure followed here is always based 
on taking appropriate limits ( t -*• 0, e ■* 0, etc. ), and whether the resulting 
equations describe sharply distinguished physical regimes encountered in 
practice must depend on detailed numerical calculations and experiments. 

II. FORMULATION OF PROBLEMS 

Consider a two-blade rotor advancing at a constant velocity V and 
rotating at angular velocity perpendicular to the direction of advance. 
For our purposes, the number of blades is not important. Attention is 
restricted to the outer few chords of the advancing blade. In the near 
flow field, close to the advancing tip region, the flow is not influenced (in 
a first approximation) by other portions of the blade system. Boimdary 
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conditions on the blade are formulated for arbitrary planform, although 
special care may be required if the planform has some unusual shape, as 
for one with a kink. We shall not consider the lifting problem because too 
little is known about rotor wake properties, this being particularly true 
for high speed rotors. It will also be assumed that the flow is inviscid 
and isentropic. A velocity potential therefore exists for the flow described 
in an inertial coordinate system. We take this system to be one at rest 
relative to the vmdisturbed air. 

The potential equation for (p is an inertial reference system is 

<?tt + (W)^l = (p (1) 


where c is the local adiabatic sotand speed, 
c and is, 


n 


i (V^)" -t fri 



Bernoxilli's equation, relating 


(2) 


where c is the sound speed in the undisturbed gas. 

00 

In a fixed Cartesian (x, y, z) coordinate system, let the blade geo- 
metry and location be described by F(x, y, z, t) = 0. The condition of flow 
tangency at the blade surface is 


+ V<p- VF = 0, (3) 

and we require that (P^ and V<P go to zero at infinity. 

In all cases of interest here (including nonlinear problems), the 
third term on the left hand side of equation (1), the cubic term, is negli- , 
gible. That term will be neglected in all further discussion. 
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It is convenient to reformiilate the problem in a coordinate system 
attached to the blade, with the origin at the center of rotation. The velo- 
city of advance V will be taken in the direction of the negative x-axis of 
the fixed reference frame. Coordinate z for both fixed and moving frames 
is normal to the plane of rotation, and the angular velocity of rotation n 
points in the positive z-direction. In the rotating frame, the blade -attached 
x-axis is in the chordwise direction and the y-axis is along the blade span. 

We now transform equations (l)-(3) and rewrite them in a coordinate 

■V. 

system attached to the blade. For a moment it will be necessary to distin- 
gtiish between time and space coordinates of the inertial frame and time and 
space coordinates of the moving frame. Let primed variables refer to the 
inertial system, and unprimed variables refer to the blade -attached system. 

I III _ 

Let r = (x , y , z ) and r = (x, y, z) be position vectors in these two sys- 
tems. The two reference systems are related by 

r' = V t + R(t ) 7, 

I 

t = t, 

where R(t ) is the rotation matrix 

cos nt - sin nt 0 

R(t) = sin f2t cos Sit 0 

u 0 0 1 

The inverse transformation is 

— _ii_» _ii_i 

r = R (t ) r - R - (t ) V t 

I 

t = t . 
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where R ^ is the inverse of R. For the purposes of equations (l)-(3), the 
gradient operation in the inertial system may be replaced by the gradient 
operation in the moving system. Time differentiation in the inertial system 
transforms as 


I 

at 


^ ^ a^ 

at ax. 


with a sum over the i- subscript. It is easily shown that 


(4) 


= - [R(0) r] - [R~ \t) V] = - (fl X r + V) 

at ^ 


where R(0) is dR/ dt evaluated at t = 0. In a more convenient notation, 
equation (4) becomes 


dcp _ dcp 

at 


^ - (jTx7+ V) . V(P 


at 


(5) 


The velocity V in equation (5) must be expressed as it appears to 

I 

an observer in the moving frame. In the fixed frame, V = - i V. In the 
moving frame. 


V = - i V cos n t + j V sin t. 

An observer in the moving frame sees a free- stream velocity -(O x r + V). 
Let 


a = - (n X r + V) . 


Then equation (5) is 


d<P 

I 

at 


fr- 


V(p 
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To traJisform equations (l)-(3) to the moving frame, we first note that a 
second application of equation (5) gives 

<P = <p + a • V • a+2a* V<^-Qxa* V<P - V • V<P . 

I , tt t 

^ t 

The second term in equation (1) is 

-V (V = (-^ + I • V) (V(/>)2 
at 

= 2 W • + 2 a . V V(|P • v<p 

2 2 2 

The right hand side of equation ( 1), c V requires c from Bernoulli's 
2 

equation (2), and c must be expressed in the moving frame. Bernoulli's 
equation transforms to 

c^. = c^ - (7 - 1) + I • V(p + I (V(p)^] . 

The complete equation for (p written in the moving frame is now 
+ 1 . V V(p . a + 2 I . Vcp^ - n X I • V(p -V . V(p 
+ 2 . V (;P^ + 2 a • V .V<p • 

= { c^ - (7 - 1) [</>. + I • V(/? + -i (V(pfy } (6) 

00 I ^ 

Equation (6) is the exact potential equation written in a coordinate system 
that translates at constant velocity V and that rotates with constant angu- 
lar velocity S2 . As far as equation (6) is concerned, V and S2 are arbitrary 
vectors as they appear in the definition a = - Jlxr-V. Later, V and 
will be specified as vectors appropriate for a rotating blade problem. 
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It is very usefvil to simplify equation (6) by non-dimensionalizing 
and scaling it in a way that is appropriate for the particiolar problems at 
hand. This simplification requires a choice of suitable length and time 
scales, and a further specification of the relative magnitudes of all non- 
dimensional parameters that enter the problem. The choice of length and 
time scales will in turn depend on what part of the flow field one is inter- 
ested in. For the rotor problem, there are several distinguished flow 
regions. Some of their features are listed here. 

Region 1. This is the region near the blade and sufficiently far inboard 
from the tip that tip effects are not important. It is the strip theory or 
blade element region. If the blade tip Mach number is subsonic or tran- 
sonic, equation (6) in Region 1 will become linear and quasi- steady* that 
is, all nonlinear terms and time derivatives may be neglected. In addition, 
the flow in this region becomes approximately two-dimensional. Span 
coordinate and time appear only as parameters in the problem. The length 
scale is blade chord and time scale is \ 

Region 2. This is the region near the blade surface and near the blade tip. 
If the Made tip Mach number is subsonic (not transonic), equation (6) will 
be linear, quasi- steady and three-dimensional. The length scale is dis- 
tance in chords measured from the tip and time scale is \ In a first 
approximation, the blade will appear to be semi-infinite as viewed from 
the tip. When a quantity (e. g. , <P ) is calculated in the tip region, then, as 
the calculation marches inboard from the tip, the results must match those 
obtained in Region 1. 

If the tip Mach number is transonic, analysis of equation (6) in the 
tip region is more complicated. The equation will be nonlinear- in spatial 
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derivatives, and may be quasi-steady or essentially unsteady, depending 
on certain limiting cases to be discussed later. 

Region 3. This region will be called the first far field region. The char- 
acteristic length scale is blade span. For subsonic flow, an observer in 
this region would view the blade as a line rotating about its mid-point. 
Disturbances measured in this region would be considered as caused by 
point singularities continuously distributed along the blade line. Equation 
_ (6) becomes linear, three-dimensional and unsteady for both subsonic and 
transonic Mach niimbers. The inner boimdary of this region is the com- 
bined flow of Regions 1 and 2; the outer boundary is the field at infinity. 

2 2 2 

Region 4. This is the distant far field, defined by letting x + y + z oo . 
It is the region appropriate for far field acoustic calculations. At any point 
in this far field, the local flow is described by retaining all linear terms in 
equation (1). Time scale is ^ and distance is a characteristic wavelength 

We are primarily interested in Region 2, the tip region. A prelim- 
inary simplification of e'quation (6) will now be made, and it can be verified 
later that terms dropped here are negligible. All first order derivatives on 
the left-hand side of equation (6) may be neglected. The term 2 
may be neglected compared to the term 2 a • because a is essentially 

a first order term and V<p is the perturbation velocity. On the right hand 

12 2 
side, the terms and — (V^) in the square bracket that multiplies V may 

be neglected. With these observations, the basic equation (6) becomes 

+ a . V V(p • a + 2 a • + 2 a • V '^<P ' V<P 

= [c^ - (t - 1) I - VV (7) 

00 
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The condition that the flow be tangent to the blade surface is easily- 
expressed in the moving coordinate system, where the blade surface ap- 
pears stationary. Let blade geometry be defined by z = g(x, y). Then the 
blade boundary condition is 

<P^= (- Q x7 - V + W) * Vg,. 

evaluated on z = g(x, y). The term is negligible, and we use the mean 
surface approximation to get 

(x.y, 0 , t) = - (n X 7 + V) . Vg. ( 8 ) 

Scaling of equations (7) and ( 8 ) requires that they be written out in 
scalar form. Using the definition of a, and rearranging some terms, 
equation (7) expands into the following lengthy form: 

+ 2{J7y + V cos nt-) (p^ - 2(nx + V sin 

= [c^ - (fiy + V COS nt)^ - (7+1) (fiy + V cos S2t) (P^ 

- {7 - 1) (S2x + V sin nt) (p^] 

+ [2(S2x -t V sin fit) (fly -f V cos J2t) + 2(S2x 4- V sin J2t) <P^ 

- 2 (fly + V cos ^t)(P^] <P^ 

+ [c^ - (fix + V sin - {7 - 1) (J2x + V cos Ot) 

+ (7 + 1) (fix -f V sin Sit) (P \ ^ 

' yJ yy 

- 2 (n + V sin Sit) <P + 2(S1 + V sin Sit) <P 9 

y z xz ' X z xz 

+ [c^- ( 7 -l)(n +VcosSlt)<? +( 7 -l)(Sl + VsinSlt)]<P (9) 

00 y ZZ 
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The expanded form of the boundary condition, equation (8), is 

(P (x, y, 0, t) = (f2 y + V cos Q t) g - X + V sin t) g ( 10) 

z X y 

Equation (9) and (10), with a suitable boundary condition at infinity, des- 
cribe the aerodynamics of the nonlifting rotor blade. These equations will 
be further simplified by scaling them. This scaling will be done in the next 
section where the two most important problems are considered, the blade 
with subsonic advancing tip Mach number, and then the blade with transonic 
Mach number on the advancing tip. 

III. SCALED EQUATIONS 
1. The subsonic problem 

We first consider the problem of a blade with advancing tip Mach 
number that is subsonic. 

The introduction of new symbols for nondimensional scaled variables 
can be avoided by a convenient notation. The combination fit in equation (9) 
is simply replaced by t, with the understanding that from now on t stands 
forSlt. Similarly, let ^ be the blade chord, and use the substitution 
x/ j? -► x, where the x on the left is dimensional and on the right non- 
dimensional. The complete set of time and spatial substitutions is then 


nt ->■ t 

x/ j? X 

/ 

y/R - y 
z/ ^ z. 

A reference Mach number will be taken as 

M = J2R/c 


( 11 ) 


( 12 ) 


11 



and advance ratio 


U = V/nR . (13) 

The blade has reference thickness ratio t. Using the substitution schenoe 
(11), the scaled potential for the subsonic problem is 


<P 

Tfi R jf 


<P. 


(14) 


The blade has aspect ratio AR = R/ i. Define the parameter e 


e = AR‘^ = i/R , 


(15) 


which is generally a small number for a rotor blade. We now use the substi- 
tutions ( 1 1) - (15) in the differential equation (9). All nonlinear terms may be 
neglected in the subsonic tip problem. The scaled differential equation is 
therefore 

M^e ^ + 2M^e (y + cos t) - 2M^e ^(e x + p sin t) (p^^ 


= [ 1-M^(y + u cos ^ 2M^e (e x + p sin t) (y + p cos t)<P^ 


+ 6 i-M^(e X + p sin t)^]^ + (P 

'■ ^ yy zz 


(16) 


A properly posed problem for the tip region requires some y-derivatives in 
the differential equation. The scale length for this region is distance s in 
chords measured from the tip toward the center of rotation; 


1-y 
s = — ^ 
e 


y = 1 - e s 


(17) 


Equation (16) becomes 
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M^e + 2M^e ( 1-e s + p cos + 2M^e (e x + ju sin 

= [ 1-M^( 1-e s + y cos t)^]</’ - 2M^(e x + p sin t) ( 1-e s + ^ cos t)(P 

XX xs 

' + [1-M^(ex + M sin t)^]<P , + (P. (18) 

s s z z 

Keep M, p, x, s, z and t fixed. Let e -► 0. Equation (18) becomes 

[1-M^(1 + p cos t)^]<P - 2M^ p sin t (1 + ju cos t)(P 

XX xs 

+ (1-M^Ai^ sin^t)</> + (P = 0 (19) 

s s z z 

In the boundary condition, equation (10), let blade geometry be fur- 
ther specified (in dimensional form) by 

z = g(x, y) = i Tf(^, ^) 

where t is a reference thickness ratio. Nondimensional scaled geometry is 

z = T f(x, y). 

Equation (10) transforms to 

(x, y, 0, t) = (y + /Lt cos t)f - e (e X + ju sin t)f . (20) 

Z X y 

Replacing y by 1-e s and letting e -> 0 gives 

' (x, 1, 0, t) = ( 1 + /u cos t) f (x, 1) + sin t f (x, 1) (21) 

The differential equation (19) and boundary condition (21) describe the tip 
flow for all values of time, that is, for all azimuthal angles around the rotor 

disk. The scaled problem is also quasi- steady, with time appearing only as 

' 

a parameter. Note also that the problem posed by equations (19) and (21) is 
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exactly the same as the problem of the flow in the tip region of a fixed yawed 

I 

wing in subsonic motion, moving with forward velocity V = -i (12R + V cos Ot) 

I 

+ j Vsinflt relative to the \andisturbed air. 

It has been assumed so far that the advance ratio p is not small. If 
p were small, say of the order of thickness ratio t, then it would be incon- 
sistent to retain terms of order p and smaller in equations (19) and (21). If 
p = 0 (t), then all terms involving p must be neglected, and time no longer 
enters the problem at alll for small p , the flow in the tip region of an ad- 
vancing rotor with subsonic tip Mach number becomes equivalent to the flow 

I 

in tip region of an unyawed fixed wing moving at constant velocity V=-i V 
relative to undisturbed air. 

While equations (19) and (21) describe the tip flow for the subsonic 
rotor, they are not convenient for use in numerical calculations. A single 
set of equations that describes the flow in the tip region, and which remains 
valid inboard from the tip through the strip theory region, is found by in- 
spection of equations (19) and (21) to be 

[1-M^(y + p cos t)^l<i^’ + 3M^e p sin t(l + p cos t)^ 

^ ' .1 jQj 

+ e%l - M^p^ sin^ t]<? + = 0 ' (22) 

^ ^ yy z z ' 

<P (x, y, 0, t) = (y + /u cos t)f - e p, sin t f , 
z X y 

2. The subsonic far field 

In numerical calculations, some condition must be imposed on a spa- 
tial boundary that is, in some sense, far from the blade. There are two 
possibilities. The first is to set (p and its derivatives equal to zero on 
a boundary that is some number of blade chords distant from the blade itself. 
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The second alternative is to derive an asymptotic far field representation of 
the solution, and impose this representation on a boundary that is again suit- 
ably far from the blade. This second approach has been used by Murman and 
Cole (2) for a two-dimensional airfoil in transonic potential flow, and extended 
by Klunker (3) to the case of a three-dimensional fixed wing in transonic flow. 
For the subsonic rotor problem, the asymptotic far field is most easily de- 
rived in an inertial coordinate system, where the potential satisfies 

= c V (P . The far field is then represented by the radiative part of a 
potential field that is caused by a system of retarded doublets continuously 
distributed along the span of the wing. The strength of the doublets will vary 
along the span, and also with time for an advancing rotor. When this repre- 
sentation has been found in the inertial system, it can be transformed to the 
moving frame and used as an outer boxmdary condition for equations (19) and 
(21). It is not clear, however, whether the additional labor involved in find- 
ing and using an accurate far field representation is worth the effort. Numeri- 
cal experiments with two- and three-dimensional fixed wings indicate that 
simply setting <p =. 0 on the outer numerical boundary leads to sufficiently 
accurate results. 

3. The transonic problem 

The starting point for the nonlinear transonic problem is equation (7), 
and its expanded form (9). We are again guided by the close analogy between 
the flow near the tip of the rotor and the corresponding flow near the tip of a 
fixed wing. The preliminary scaling of spatial coordinates is exactly the 
same as for the fixed wing problem. Thus, in equation (9) we make substi- 
tutions 
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n t -* t 

x/ i -►X 

y/R — y 

1/3 /* 

T 7.1 i -* Z 

<P/j?^\b.s. -* q> 


After deleting some terms by inspection from equation (9), we get 


M e 


.2 e 


2M^e^ 


2 ^ + 2 M 2^3 (y + W cos t)<p^ 2/3 ^ ^ ®1” 1)*^. 


yi 


_ [ -1-M (y^i^cos t) _ M^(y + p cos t)<P^ <P^ 


2M.^e 


+ 2/3 ' t) (y + P cos t)<P^y ■'■^73 ^yy + 


zz 


(23) 


The blade boxmdary condition is 


^^(x, y, 0, t) = (y + (u cos t)f^ - e(ex + /u sin t)fy. (24) 

In contrast to the subsonic problem, some as sumption. must now be 
made about the relative orders of magnitude of e = AR ^ and t. A typical 
blade has thickness ratio of the order . 06 to 0. 1 near the tip, and aspect 
ratio of order 20 to (in an extreme case) 10, Therefore we assume the 
limits T -*0, e -* 0 and^/ e = 0(1). We first assume that the advance ratio 
is of order 0. 3 to 0. 4 and is held fixed when any limits are taken. 

When the flow is transonic only on the advancing side of the rotor 
disk, there is some finite azimuthal angle measured about t (or 12 1) = 0 
within which the flow is transonic and outside of which the flow should be- 
come subsonic. It is therefore necessary to scale time in a suitable way. 
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Also, when the flow is transonic there is some characteristic region, mea- 
sured in chords, near the tip when transonic conditions exist and inboard of 
which the flow becomes subsonic and follows strip theory. It is therefore 
also necessary to appropriately scale the span coordinate y to obtain a pro- 
perly posed boiondary value problem in the tip region. Following Reference 
(1), we introduce a new span variable: 




(25) 


Introduce a time scale parameter 6 by 


t = 5 T, 


(26) 


with 6 to be determined. Use of equations (25) and (26) in (23) gives 

2 2 2 
M e , 2M e „ e ^ 

■2737^ ^TT + ~2737 (1 - -173 S + M cos 6T)<? 

T ' 0 T ' 6 T 


= [ 


+ 2 (^^ + M 

S + 

2/3 


1 - M^(l— S + p cos 6T)^ 


(y+l)M^(l--^ S+^ cos 

T ' 


2M 


Hy ' g (c^x + (ji sin 6T) (1 - ^ S + p cos 6T) 


+ + ^ . 

So zz 


(27) 


The boundary condition is 


1 / 3 / 


<^^(x, 1 - S, 0, 6T) = (1 - -j y ^ S + p cos 6T)ljj + t' (ex+psin6t)f 


The maximum Mach number of the free stream relative to the blade occurs 
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when T = 0, and is . W = M(1 + ^). As for a fixed wing, we assume 

(1 - TJl^) T = 0(1). Consider first the coefficient of in equation (27), 

and in particular consider the part 



Let S and T be 0(1), let e -*0, t -* 0, €/ t = 0(1), (1 - ^ = 0(1), and 

for the moment assume also 6 -*-0 as t-*0. Equation (28) becomes 



The second part of the square bracket that multiplies becomes, in the ' 

2 

above limit and with x = 0(1), (y +1) M (1 + which provides the usual - 

nonlinearity. We assume this term is 0(1). 

The magnitude of the last term in equation (29) is determined by the 
-l/ 3 

limiting value of 6 t ' • , with both 6, t— ►O. Examine the coefficient of 
in equation (27) in the above limits. The coefficient of becomes 

- 2M^ u 6 (1 + <i) T. 

-l/ 3 

If 6 T ' becomes unbounded as t— ►0, the (P ^ ^ terms would dominate 

XO XX 
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all other terms in equation (27), and the scaling procedure followed so far 

-l/ 3 

would be incorrect. Alternatively, if 6 T -*■ 0 as t 0, inspection of 

equation (27) shows that the term alone dominates, leading to another 

-l/ 3 

useless result. The conclusion is that 6 should be chosen so that 6 t ' =0(1) 

as T -*■ 0, and we therefore take 


6 = 



Combining these results,, the asymptotic form of equation (27), valid in the 
tip region on the advancing side, is 

2 

2M^(1 + ^)^ + 2M^(1 + u)^s+ M^(l + p) p 

T 

- (y + 1) M^d + (U) < P ^ < P ^ 

- 2M^ M (1 + p) T + <;pgg + (30) 


while the blade boundary condition is 

<P^(x, 1, 0, 0) = (1+ p) f^. (3 1) 

The original choice of the relative order of magnitude of € and t 
was based on the observation that e = 0 (t) for a typical rotor blade. How- 
ever, it now emerges from equation (27) that, analytically, a necessary con- 
dition for the flow to be transonic and nonlinear is that precisely this relation 
between e and t shoizld hold. 

The transonic problem defined by equations (30) and (3 1) differs from 
the subsonic problem in two essential respects: the transonic problem is 

nonlinear (in the usual way, and as is to be expected), and it is essentially 
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an unsteady problem with the coefficient of = 0(1). Thus the simplifying 
feature of quasi- steady flow does not appear. This latter result naturally 
complicates numerical calcvilations because solutions must be obtained in a 
four-dimensional (x, y, z, t) space. The term in however, again disap- 

pears. The hypersurface T = constant (or t = constant) is therefore charac- 
teristic (the necessary and sufficient condition for T = constant to be charac- 
teristic is that the term in not appear in the differential equation). 

Notice that the tip flow is transonic, unsteady and nonlinear over a 

1/3 

sector of the rotor disk of order t (or Qt) ~ t' , with t measured about t = 0. 
If T= 0. 1, the azimuthal extent of the transonic region is about 54°, or 27° 
either sice of t = 0, 

While equations (30) and (3 1) are a first approximation to the flow in 
the tip region, they are not convenient for numerical calculations. Following 
the same procedure that was used for the subsonic problem, we retain only 
the appropriate terms in equations (23) and (24) and find 

2 2 

2M^ (1 + ju cos ^ (y+1)m2(1+|u cos t)(p^YP^ 

T T 

2 2 

+ Ml + M cos ^ 

^ = (y + W cos t)f on z = 0. 

Z X 

These two equations describe the flow in the transonic tip region and remain 
valid inboard through the subsonic strip theory region. 

4. Reduced frequency 

A fixed wing with chord I moving at constant velocity V and oscil- 
lating with cingular frequency has associated with it a nondimensional re- 
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duced frequency lo = / V. This parameter is a measure of the number of 

oscillations the wing goes through as a fluid particle goes from leading to 
trailing edge in its mean motion. Unsteady flow is generally characterized 
by large u and quasi- steady flow by small u. Usually u is multiplied by 
some other parameter or parameters in the governing differential equation 
. and boundary condition, and appropriate limits must be taken to specify what 
is meant by large or small u. In the rotor problem, a reduced frequency in 
precisely the form fii/ V does not appear, but a similar parameter can be 
interpreted as a reduced frequency. The reference velocity near the tip on 
the advancing side is JIR + V. The reference transit time of a particle in its 
motion from leading edge to trailing edge is i. (J2R + V) The nximber of 
radians traversed by the tip in this time is 

“ = S2R + V • ” 


i/R 

1 + H 


This parameter may be defined as the reduced frequency for an advancing 
rotor. Since (x is never large, the reduced frequency is effectively the in- 
verse of the blade aspect ratio e. Reference to the basic differential equa- 
tion (23) indicates that e does indeed play the role of reduced frequency 
wherever it multiplies time derivatives. 


5. The transonic far field 

The specification of a far field boundary condition in the transonic 
problem is considerably more complicated than for the subsonic problem. 
This is particularly true if an asymptotic representation more accurate than 
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<P = 0 is to be prescribed on a boundary some chords distant from the blade, 

2 2 2 

If a distant far field solution based on the limit x +y +z — ►oois reqmred, 
then finding such a solution is equivalent to finding the radiative part of the 
solution at infinity, and getting this solution means solving the acoustic prob- 
lem for the rotor. One method for analyzing the acoustic far field for a body 
in transonic motion has been described by Ffowcs Williams and Hawkings (4), 
but it has not yet been applied to rotor problems. The method is rather com- 
plicated and involves the near field transonic solution in an integral repre- 
sentation of the far field. The essential difficulty here is that an element of 
blade surface in transonic motion causes a disturbance that is not compact, 
and the far field of the disturbance cannot be represented as a simple point 
multipole, as it could be for subsonic motion. These remarks indicate that 
without a much deeper analysis, there is no alternative, at the present time, 
to using the condition <P = 0 at a far field boundary in numerical calculations. 

6. The transonic problem with small advance ratio 

It has already been noted that for a subsonic blade at moderate ad- 
vance ratio the flow is quasi- steady and linear; and when advance ratio is 
small, time does not enter the problem at all; if, say, p = 0(t), then the 
effect of forward motion becomes a second order effect, and the first ap- 
proximation to the flow is simply the flow over a hovering blade. 

We now consider the transonic problem for small advance ratio and 
determine the appropriate equations in the limit p 0. 

Refer again to the basic differential equation (23) and consider the 
first term in the square bracket that multiplies 
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1 - M^(y+u cos t)^ 
2/3 

T 


1-M^{1 - S + /u cos 6T)^ 

271 


E 2M^(l+/u) ^ S 

2 2 2 
+ (1+M - S)m(1-cos 6T) - (1-cos 6T)2 


jn 


S. 


(32) 


2/3 

We first consider the case for which ju — ► 0, /n = 0(t ). Set the 

time scale factor 6 = 1 in equation (32), and revert to 6T = t. Let S = 0(1), 
t= 0(1), e/r= 0(1), lu = 0(t^^^), e — ► 0, t — ► 0, n — 0. The asymptotic 
form of equation (32) is 

l-M^(t+w cos t)^ 1-M^ „,,2 u 

^/f • jrr - 2M ^ + 2M - S 


+ 2M^ ^ (1-cos t) 


Time derivatives become negligible, and the asymptotic forms of the differ- 
ential equation (23) and boundary condition (24) are 


- 2m 2 ^ t 2 m2 1 S + 2U^-^ (1-co. t) - (y*l)M\]<P^ 


+ ‘^SS + = O- 


^ (x, 1, 0, t) = f . 
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The cross derivative is also negligible in this limit, its coefficient being 

1/3 ^ 

0(t ). The flow is therefore quasi- steady and sweep-back effects are of 

second order. The flow is, however, transonic aind nonlinear around the 
entire rotor disk. 

The differential eqviation and boundary condition in forms suitable 
for computations from the tip inboard through the strip theory region are 






(p +(p = 0, 

yy zz 


(33) 


(P^ = y f^, on z = 0. (34) 

2/3 

There is also a limiting case based on the limit p / t ' 0. It is 

apparent from equations (33) and (34) that in this case time derivatives, 
sweep back effects and time itself do not enter the problem in both the tip 
region and the strip theory region. In this limit, the flow over the rotor 
from the tip inboard is therefore equivalent, in the first approximation, to 
the flow over a hovering rotor with tip Mach munber M. 

7. The very lar ge aspect-ratio blade 

A final limiting case emerges naturally from equations derived in 
previous sections. It is based on the limit t 0, e -♦ 0, e/r-^O and corres- 
ponds to a very large (and perhaps impracticably large) aspect ratio blade. In 
this limit, time derivatives and sweep-back effects disappear, and also the 
effect of a mean Mach number gradient along the blade span becomes negli- 
gible. For a hovering blade, this limit implies that the rotor problem as 
described in a blade fixed coordinate system is entirely equivalent to the 
problem of a fixed wing with free- stream Mach number J2R/ c. For the 
advancing blade, there are several flow regimes associated with e/r-» 0, 
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depending on the order of magnitude of the advance ratio p. In all cases, 
the flow over the very large aspect ratio blade is analogous to some fixed 

f 

wing problena. Further details are imnecessary here, since each limiting 
case can be read directly from the differential equation and boundary condi- 
tion corresponding to a particular choice of /u. 


IV. PRESSURE COEFFICIENT 


The pressure coefficient C is based on the local chordwise free 
^ P 

stream velocity at each span station. In nondimensional form, this velocity 
is y + p cos t. The leading terms in the expression for C^ are 

P - P^ 


C 

— £ =1 2 2 2 
T -y p n R (y+p cos t) T 

L* 00 


2 


y + p cos t 


2e p sin t 
+ ^ 


(y+p cos t) 


2 • 


(35) 


For subsonic flow in the tip region this becomes 


C Z(P 

X 

T 1 + P COS t 


2 p sin t 


(1 + p cos t) 


2 • 


(36) 


Strictly speaking, if p = 0(1), the second term involving the spanwise de- 
rivative ^ must be retained. While this term may be neglected in a sufficiently 
small sector about t = 0, it should be noted that there is no characteristic azi- 
muthal angle on the advancing^ side for a subsonic tip (in contrast to the advanc- 
ing transonic tip, for which a characteristic azimuthal angle about t = 0 does 
play an essential role). The terms in in equation (35) and ^ ^ in equation 
(36) should therefore be retained in a subsonic calculation when p = 0(1). 

The second term in equation (3 5) is never important in transonic 
flow. Therefore, for this case 
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( 37 ) 


C 2(P 

V _ _ X 

2/ 3 " y + U cos t 


Equations (35) and (37) represent from the tip inboard through 
the strip theory region for subsonic and transonic tip Mach numbers. 


V. CONCLUSIONS 

The most difficult problem for numerical calculations occurs when the 
advancing tip in transonic and advance ratio is about 0. 4. Under these condi- 
tions, a shock wave begins to form as the blade enters the advancing side of 
the rotor disk. It grows to some maximum strength in the advancing region 
and then decays as the blade retreats. Both the shock wave and blade thickness 
cause weak waves to be propagated very slowly, relative to the blade, into the 
upstream region. The slow propagation speed of these waves causes the essen- 
tial unsteadiness in the flow and complicates the governing equations. If the 
tip is transonic and advance ratio is 0, 4, a blade of radius 20 feet would have 
a forward velocity of about 190 knots and 375 rpm at sea level conditions. 

These conditions could be encoiintered in cruise or maneuver for a high per- 
formance helicopter. 

The computational problem is much simpler when advance ratio is small. 
The flow becomes quasi- steady, sweep-back effects are negligible and the flow 
is nonlinear and transonic around the entire rotor disk. Computer programs 
for a hovering transonic rotor can be used here because time enters only as 
a parameter. 

If advance ratio is very small (of the order of thickness ratio or smaller), 
sweep-back effects, time itself and therefore time derivatives do not enter the 
problem in both the tips and strip theory regions. In a first approximation, this 
problem is equivalent to a hovering rotor. It differs from a fixed wing problem 
only by the mean Mach number gradient along the span. 
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other limiting flow regimes for subsonic and transonic tip Mach num- 
bers are discussed in the main text. These regimes have interesting physical 
features, but they are of secondary importance in their numerical, aerodynamic 
and acoustic properties relative to the high speed, high advance ratio blade. 
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